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The solution to an initial value problem for the flow of a buoyant inviscid incompressible fluid 
in three-space is given in terms of a fundamental solution. The initial values and the distribution of 
source strength can be quite general subject only to mild restrictions on smoothness. The fundamental 
solution and its associated velocity and displacement fields are given explicitly and in some detail. 
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Introduction 

This article presents a solution to an initial value problem in the theory of stratified fluids. 
Thus it is a contribution to the theory of buoyant flow and internal gravity waves and also to the 
closely related topic of flows in rotating fluids. 1 A system of two linear second order equations in 
three space dimensions and time govern the problem. The equations describe the flow of an in- 
viscid incompressible fluid in which an arbitrary distribution of source strength is prescribed. 
The approximations that the density and the buoyancy frequency are constant are used. A funda- 
mental solution of the system is used to solve the initial value problem in terms of the source 
distribution and the initial conditions. The fundamental solution is given explicitly; and a number 
of its properties are noted, including its velocity and displacement fields. 

Section 1 gives a summary sketch of the system of equations governing the flow. They result 
from a first order perturbation from static equilibrium of the conservation equations for momentum, 
density, and mass. Needless to say, these equations agree in all respects with those given in 
standard works such as [8]. 2 Two technical differences are worth noting. Kirst, the unknown 
functions in the system with which we work differ somewhat from those most often seen; and 
second, we prefer to deal with a pair of second order differential equations rather than combining 
them into a single fourth order equation. This section concludes with the statement of an appro- 
priate set of initial conditions for the problem. 

We start section 2 by writing down the full initial-value problem in the exact form we need. 
Then in section 2.1 we give the fundamental solution and an interpretation of it as an impulsive 
injection of fluid; the corresponding velocity field is also presented. A verification of the solution 
is sketched. In section 2.2 the solution of the problem with a prescribed source distribution but zero 
initial conditions is given. Then the initial conditions are accounted for, and everything is put 
together. Some of the properties of the solutions are verified in this section, but other technical 
details will appear elsewhere [3]. 

Properties of the solutions are examined in section 3. The structure 1 of the fundamental solution 
shows a natural decomposition of the flow into two parts, one the flow of nonbuoyant ideal fluid of 
constant density and the other a modification that gives the influence of buoyancy. The first part 
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responds immediately to the injection, the second is slow and long lasting. We give especially 
simple expressions for the fundamental solution directly above and below the singular point and 
in the horizontal plane with it. 

The solutions of the initial-value problem are then examined in light of the properties of the 
fundamental solution. We note that, in addition to an oscillating source-like effect, the initial condi- 
tions produce a term that rings without changing spatial distribution. Under certain conditions this 
can be the total effect of the initial conditions. In other circumstances, when the initial conditions 
are confined to a bounded region, we show that the influence of the initial conditions dies off at 
large distances faster than that of a source distribution with nonzero total strength. 

The last section discusses very quickly the displacements that result from the fundamental 
solution. These are purely radial and are found from an exact integration of the velocities. 

A summary of portions of this work has been presented in [2]. 

1. The Equations of Motion 

We start from a usual system of equations for an inviscid incompressible fluid in a constant 
gravity field. These state the conservation of momentum, density (incompressibility), and mass: 

pDulDt + gradp + gpk = 0, (1.1) 

Dp/Dt = 0, (1.2) 

p div u = /x, (1.3) 

where k is the vertical unit vector (in the ^-direction), and p is a prescribed distribution of source 
strength. The rest of the notation is standard. These equations are to be linearized about a state 
of stable static equilibrium described by u=Uo = 0, p = po(x^), and p = Po = — g f*' po(z)dz. For 
stability p () is required to be positive and monotone decreasing with x%. The perturbed flow is taken 
in the form u= Uo + u, p = po + p, and p = p {) -\-p. These are put in (1.1), (1.2), and (1.3), and quan- 
tities of order higher than the first are discarded to obtain the linear equations 

podu/dJ + grad p-f gpk = 0, (1.4) 

dp/dt + u 3 dpoldx:]= 0, (1.5) 

p div u^p. (1.6) 

Now we drop the 's, and henceforth all the quantities will be the perturbations except for p , 
which designates the prescribed equilibrium density. From the horizontal components of (1.4) it 
follows that the vertical component of vorticity does not change with time so that for a variety of 
flows, including all motions that start from rest, the horizontal velocity components derive from a 
potential. We assume this to be the case. For convenience we write this in the form 

1 <v 

Ua= — -r* 1 , a- 1,2. (1.7) 

Pi) OX a 

Also from (1.4) and (1.7) it follows that p can be expressed by 

-p=d<pldt, (1.8) 

and that the equation of motion in the vertical direction takes the form 
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at \ 0X3 

The equations (1.9) and (1.5) imply 



tpou 3 -^-)+gp = 0. (1.9) 



and (1.6) takes the form 



d 2 / 1 d<p\ g dp 

Tr* "3 — t~ — ~y~ "3=0. 

d£ 2 \ p () d# 3 / po a^s 






Now we introduce the notation n and w*. The buoyancy frequency n is defined by ra 2 = — (glpo)dpol 
dxz, and w* is given by 



_J_i£. 



m;* = m 3 f- (1.10) 

po 0*3 



We write the preceding pair of equations in the more convenient forms 



and 



a*} 



+ n* „,*+— -H = 0, (1.11) 

ot* J po 0x3 



+TT+P0T— — 7^ ) + Po— = /*. 1.12 



Note that in general n 2 and p<> depend on jc. { . 

The coupled system (1.11) and (1.12) governs the quantities </? and w*, from which in turn the 
eqs (1.7), (1.8), (1.9) and the definition (1.10) determine the velocities, density, and pressure. How- 
ever the system must be completed by appropriate boundary conditions for <p and w* on the pre- 
scribed boundaries, a surface condition on free surfaces, and by initial conditions. 

Here we shall be concerned with the situation in which all of the boundaries are far away from 
the region of interest and their influence can be neglected. In this case only initial conditions remain 
to be considered. These fall on w* alone since in (1.11) and (1.12) w* is the only quantity to be 
differentiated with respect to t. They can be prescribed by 

w*(*,0+ )=/(*) and — -(x,0 + )=g(x) (1.13) 

at 

where /and g are given functions. The second of (1.13) is equivalent to a statement of the initial 
perturbation density as can be seen from (1.9) and (1.10). 

It may seem strange at first to give initial conditions on just w* and dw*ldt, but this is precisely 
what the system (1.11) and (1.12) requires. This is to be compared with a similar circumstance in 
the linear theory of surface gravity waves on an ideal fluid in which the initial values of the potential 
and its time derivative on the free surface are all that is required. 3 In each case the reason is much 

1 See |7| ( lhapter 6 lor a clear exposition. 



Q1 

211-694 O - 76 - 6 OA 



the same: an equation of motion, in our case (1.12), determines the remainder of the solution at 
t = + from what is given. 

2. An Initial Value Problem 

Here we present a solution to a general initial-value problem for a fluid in which the density 
and the buoyancy frequency are treated as constant, i.e., the flow is described by (1.11) and (1.12) 
in which p and n are taken to be constants. The region of the solution is all of three-space. 

We shift the time scale to t = nt, set g = g/n, note that ft = njn = 1, and then drop the ~'s for 
simplicity; the problem is thereby put in the form 



dt l ) p OX 



3 



V72 , dw * 



> *,*e/? 3 X(0,oo) 



(2.1) 



where V 2 is the three-dimensional Laplacian. The initial conditions become 

dw* 



w*(x, 0+ )=/(*), 



at 



(x 9 0+)=g(x). 



(2.2) 



Our approach to the solution of the problem stated by (2.1) and (2.2) is the determination and 
systematic use of a fundamental solution of (2.1). 

The importance of a fundamental solution derives from the fact that for the equations and 
region for which it is valid, it reduces the solution of inhomogeneous initial-value problems to 
integration. Further, by a variety of more-or-less standard devices it often can be modified to apply 
to a number of other regions. Equally important, an analysis of the behavior of a fundamental solu- 
tion can yield considerable information on the properties of whole classes of solutions, any of which 
might be difficult to obtain explicitly. In addition, it is a usual starting point for a formulation of 
initial-boundary-value problems in terms of integral equations. 

2.1. A Fundamental Solution 

A fundamental solution for the problem stated by (2.1) and (2.2) is a singular solution <£>, iv* 
of (2.1) that gives the effect at x, t of a disturbance at y, r. Here the disturbance can be interpreted 
as an impulsive injection at time r of a unit mass of fluid of density p ( , at the point y. The body of 
fluid is supposed to be at rest prior to the injection. Since the fundamental solution will depend 
only on x — y and t — r, we need to consider only y— and r = for the moment. 

The fluid is at rest prior to £ = 0, and consequently both <p and w* will be identically zero for 
t < 0. For t ^ our fundamental solution is given by 



^ (X ' ^ = 4ddr W+ ^' /o ^ COS 0) + L/o(t cos 6)Y 



w*(x,t) = 



-1 



47rpoM 2 I dt 



T-tJi(t cos 0) + [tJi(tcos e)Y 



7.(0 



J lit) 



(2.3) 



where \x\= {x\ + x\ + %*) 1/2 , 8(t) is the Dirac 8-measure concentrated at zero, and cos 6 — x 3 l\x\. 
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The symbol * indicates convolution, e.g., l(t)*m(t) = l{T)m(t — r)dr: and J and J x are Bessel 

Jo 

functions of the first kind. The fundamental velocities u are those corresponding by (1.7) and (1.10) 
to the fundamental solution (2.3). These formulas yield by a direct calculation. 



47rpo|*| 3 [ dt 2 



8(t)+T^tJ (t cos0) + 



T-tJo(t cos 0) 

at 



Ji(t) 



(2.4) 



In the computation it is economical to differentiate the quantity in braces in <p with respect to cos 
and then cos 6 with respect to each of the space variables. 

We give elsewhere [3] a detailed and mathematically complete demonstration of the validity 
of the fundamental solution and the representation formulas that make use of it. Here we give only 
a very brief sketch. 

Note first that <£, m;*, and u given by (2.3) and (2.4) have the correct coordinate singularities for 
a fundamental solution [1], It is quite easy to show that <p and ib* do satisfy (2.1) with /x equal to 
zero for any nonzero x. In fact, the first of (2.1) is demonstrated by showing that d 2 w*/dt 2 is just the 
negative of the third component of u given by (2.4), and consequently equal by (1.10) to 

— I w*-\ — )• Note that the d(t) appears in the second differentiation of iv* due to the jump 

\ podX;\J 

in div*/dt dtt = Q. The verification of the second equation of (2.1) is equivalent to showing that po div 
ii = for #t^ 0. This is accomplished by a routine calculation on (2.4) using the same device as in the 

computation of ii. The vanishing of po div li for i/O implies that the mass outflow (J(t) from the 

origin can be calculated from Q(t)=po J unds for any fixed closed surface S that contains the 

origin. Taking S to be the sphere | jc | = 1 we find Q(t) = 8(t). The computation is aided by using 

polar coordinates, observing that I Jo{t cos 6) s\n f) d0 = - I Jo(r)dr, and using the relation 1 

Jo t Jo 

[./o(0]* [Ji(t)lt] =Ji(t). Notice that only the leading term in (2.4) gives a nonzero contribution 
to 0(0- 

2.2. Representation of Solutions 

The fundamental solution (2.3) provides a representation for solutions of the initial value 
problem (2.1) and (2.2) with only mild restrictions on jjl,J\ and g. 

To begin we consider the case with / and g equal to zero. Then <p and m;* are given by the 
integral formulas » 

<p(x,t)=\ <p(x — y,t — T)fi(y,T)dydT, 

JO JRA 

w*(x, t)=\ w*(x — y, t — T)ix{y, r)dydr. 

Jo J R:\ 

The parallel formula for the velocity is 

u(x, t)=\ u(x — y, t — r)[jL{y, r)dydr. (2.6) 

Jo JRA 

The verification that the velocity derived from (2.5) by the formulas (1.7) and (1.10) is that 
given by (2.6) and that the equations (2.1) are satisfied by (2.5) is given in [3]. The arguments there 



(2.5) 



4 See [9] p. 380. 
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are principally of the type used in [5] to demonstrate the differentiability properties of generalized 
domain potentials for elliptic partial differential equations. A suitable set of conditions to impose 
on jut are that it be piecewise continuous in t for t ^ for all x and satisfy a uniform Lipschitz con- 
dition in x for all t ^ 0. These are sufficient for the proofs in [3], but certainly not necessary. 

On the other hand it is comparatively easy to see directly that w* given by (2.5) satisfies 
a * 
w * (x, + ) = (x, + ) = 0. In fact, the expression for w * in (2.3) shows that for small positive t 



dt 



1 



Xs 



47rp | x 



*[1 + * 2 £(cos 0, OL 



where B is uniformly bounded in t and cos 6. This implies that w*(x, t) given by (2.5) is 0(t 2 ) 
provided only that /jl(x, t) is not too rough, e.g., continuous and uniformly bounded, for small 
nonzero t. From this the required initial conditions follow immediately. 

Now suppose fi = but that the initial functions/ and g for w* and dw*/dt are nonvanishing. 
We require that df/dxz and dg/dx^ exist and satisfy the same sort of smoothness conditions given 
for fi earlier. Let h(x, t) be defined by 

h(x, t) = po[f(x) cos t + g(x) sin£]. 

Then a solution to (2.1) with (2.2) is given by 

tp(x, = - v(x — y,t — r) r— - (y, r) dy dr, 

Jo J«3 dy 3 



w*(x, t)- 



w*(x — y,t-r) " (y,T)dydr+- 



(2.7) 



dys 



Po 



and the corresponding velocity is 



u(x, t) 



T f nix- 
Jo Jr* 



y,t- 



dh h(x, t)k 

t) t— (y, t) dy dr + . 

dy.3 po 



(2.8) 



Since the integral expression in w* given by (2.7) and its derivative with respect to t vanish at t = 0, 
the function w* has the same initial conditions as h/po as required. Further, since h satisfies 
d 2 h/dt 2 -\-h = 0, the functions </> and w* satisfy the first of (2.1). The second is also satisfied with 
fi = since the contribution from the integrals in the solution is exactly cancelled by the contri- 
bution from the nonintegral term in w*. 

By linearity the sum of the solution given by (2.5) and that of (2.7) satisfy the full inhomo- 
geneous initial-value problem (2.1), (2.2); the corresponding velocities are the sum of (2.6) and (2.8). 
Thus a full solution of (2. 1), (2.2) is given by 



<p(x,t)= I I <p (x - y, t - t ) ijl (y, t ) - — (y, r ) 

■*(*, t)=\ w*(x-y, t-r) n(y, r) -— (y, r) 

Jo Jr* l dys 



dy dr, 



dy dr + 



h (x, t 



and 



u(x, t) = I u(x — y, t — 

JO J «3 



Mi(y» t) -t— (y, t) 

oy3 



c?y G?T + 



Po 



(2.9) 



(2.10) 
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3. Properties of the Solutions 

Here we examine some of the properties of the solutions of the initial-value problem (2.1), 
(2.2). We begin by pointing out some of the more evident properties of the fundamental solution 
(2.3) and the fundamental velocities (2.4). 

By inspection of the expressions for <p and u it is clear that each of them breaks naturally into 
two parts; the first is exactly that of the fundamental solution in R s for an ideal fluid of constant 
density po, and the remaining part the required modifications which account for the effects of 
buoyancy. Thus the potential and velocities can be conveniently written as <p = <p°+ if* and u = 
u°+ u*, where 



H> *(M) = ~ TZVZ\ f I: Mt cos 6) + [Jo(t cos 6)] * 



4-77 \x\ dt 



X I d 2 

u *(*,*) = A _ j ,, t \ — tj (t cos 0) + 



Ji(t) 



4<7Tp {) \x\ :i \ dt 2 



dt 



tj {) {t cos 0) 



^j 



(3.1) 



Examination of <£>* and u* shows that each of them is 0(t) at t = 0, and consequently the influence 
of buoyancy takes some time to have an important effect; however that effect is long lasting, as 
we shall see more fully later. Since the quantities in braces in (3.1) are bounded functions of cos 
for each t, it is clear that the velocity field decays in space like l/|x| 2 , i.e., just as in the ideal 
fluid case. 

As we have noted earlier, the fundamental velocity field u has vanishing divergence except 
at x = and has a net mass outflow d(t) from the origin, but the same is true of the ideal fluid part 
u°. This means that the buoyancy induced velocity u* is also divergence free and, while singular 
at the origin, has zero net mass outflow. 

The symmetry properties of the fundamental velocity u are those to be expected: axial sym- 
metry about the x* axis and odd symmetry about the origin. It is perhaps more surprising that 
u is purely radial. Since it is, and is also proportional to 1/r 2 , the corresponding fundamental 
displacements are also purely radial and can be calculated easily. 

We now turn our attention to those places where the fundamental solution and fundamental 
velocities have simple expressions. These are the x :i axis and the plane #3 = 0, i.e., 0=0, 77\ and 
0= 7r/2. In fact, we find for 0=0 and rr 



and for = tt/2 



<p-- 



-1 

477 1 JC 



6(0, w* = - 



1 



X3 



47rpo47r|#| a 



sin£, 



u 2 = 0, u 3 = 



1 



X3 



477p \x\ 



{6(0 — sin*}, 



-1 



<P = 



47t|x| 
1 



Jo T 



Utx z 



X a 



477p,) \x\ li 



6(0 



Ji(r) 



CM 
Jo — 



dr 



a=l,2; 



(3.2) 



(3.3) 



u, = 0. 



Thus, directly above and below the disturbance the vertical velocity is excited into a permanent 
oscillation at the buoyancy frequency, while in the plane the velocities it a oscillate about and tend 

1 X a 

to the values tt:' It is not really surprising that the fluid injected at t= should try to squeeze 

477po|^|' i 

into the lamina x 3 = from which it came. 
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Some of the properties of the fundamental solution that have just been noted are reflected 
immediately in the properties of the solutions (2.9) and (2.10). 

The decomposition of <p and u has its direct counterpart. In fact, we see that in general we 
may write the solution of (2.1) and (2.3) as ip = <p°+ <p*, w*, and u=u°+ u*, where 



<P°(x,t)= \ 

JR: 



-l r , , dh . 

— A ii(y,t)-—(y 

4>tt\x — y\ L dy 3 



J) } dy, 



dh 



<p *(xj) = y*(x-y,t- t) jJL(y,T) - — (y,r) \ dy dr, 

Jo Jrz L dy 3 



w*(x,t)- 



n 

JO JRs 



w*{x — y,t — r) 



p(y,v - t~ (3 / '» T ) r d y dj ^ » 

dy 3 J po 



(3.4) 



and 



47rp \x — y\ 



v>(y,t) 



dh 

dys 



,t) J dy, 



u°(x,t)= I 

JR3 

u*(^) = u*(x — y,t - t) J LL{y,r) — — (y,r) i dy dr + 



h(x,t)k 
Po 



(3.5) 



Clearly, the effect of a source distribution jjl is two-fold. First there is an ideal fluid response 
given by the expressions for <p° and u°; then there is a slower long-lasting response due to buoyancy 
given by the expressions for <£>*, w*, and u*. The latter depends, as the integrals show, on the 
history of /jl. 

Suppose that fi is a source distribution that vanishes outside a fixed bounded region G of space 
and has a total strength M(t) = J G /Ji(x,t)dx. By a standard line of argument 5 it is clear that at 
distances from G that are large compared to the diameter of G the source distribution has very 
nearly the same effect as a concentrated source of strength M located at the center of G. This leads 
naturally to an interest in a variety of concentrated source distributions, but discussion of this 
topic would lead too far afield; we shall pursue it elsewhere [4]. 

It is interesting that at large distances directly above and below G one is led by (3.2) to expect 
very small horizontal velocities and vertical velocities close to 



u 3 = 



i gg— n 



iirpo \x — y\ 



M{t)- 



■/: 



M(t) sin (t-r)dr 



where y is the center of G. The second term in the brace is exactly the response of a simple oscillator 
having the buoyancy frequency and driven by a force M{t). 

The influence of the initial conditions through the source-like term in (3.4) and (3.5) has 

dys 

the same decomposition into an instantaneous part and a slowly responding part. The terms in 

<p° and u° as well as h(x, t)lpo in w* and u* ring at the buoyancy frequency, and the latter with 

no change in shape from that given by the initial functions /and g. 

If /and g are smooth and vanish outside a bounded region G*, — dh/dy 3 has zero total source 
strength. Since h vanishes outside G*, the effect of such an initial disturbance dies off at large 
distances faster than that of a source distribution of nonzero total strength. 

It is worthwhile to note that if/ and g are independent of the vertical coordinate, then dh/dy 3 
vanishes, so the total influence of the initial conditions is given by the term h/p inw* and u>*. This 



5 See, for example, [6] (Chapter II, Article 6. 
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can be seen directly from (2.1) and (2.2) as well. In fact, it is clear that the second of (2.1) is satisfied 
with fi = by (f = and w* independent of %z. Then the first of (2.1) and (2.2) imply 

w* = f(x) cos t + g(x) s'm t= h(x, t)l po. 



4. The Fundamental Displacements 

This section is devoted to a brief examination of the displacements caused by the fundamental 
velocity field. The disturbance is taken to be at the origin, and for convenience it is scaled by the 
factor a so that the velocities are given by those of (2.4) multiplied by a. 

Since the velocity field is radial, each fluid particle is displaced radially by the flow depending 
on its original radius r before the disturbance, its polar angle 0, and time. In fact, the radial posi- 
tion r of a particle satisfies the ordinary differential equation dr/dt= u r , which takes the form 



dr 



where F is given by 



dt 47rpo r 2 



d 2 



-,F(cos 0, 0, 



(4.1) 



F(cos0, t) = 8(t)+j-,tj {) (t cos 0) + 



J- t tJv{t cos 6) 



J i(t) 



The equation (4.1) is easily integrated to obtain 

r(r o ,0;*) 



4+-^F(cos6,t) 
47rp 



i/:* 



(4.2) 



(4.3) 



in which the quantity F, defined by F(cos #, t) = \ F(cos 0, r)dr, is given explicitly by 

Jo 

F(cos0, t)=^-tj (t cos0) + I>Jo(*cos0)] * ' 

dt 



(4.4) 



From (4.4) simple expressions can be obtained for F on the x :i axis and in the plane x ;i = 0. These 
are 



F(cos 6,t) = < 



(4.5) 



COS t, 6 = 2 0,7T, 

r< J At) 

Jo(t)+t dr,2 0=7rl2 

Jo r 

47T 

For ease in interpretation we take a = — p () so that the coefficient multiplying F in (4.3) is 

one. This is equivalent to an impulsive injection at t = of just enough fluid to fill a sphere of 
unit radius. Since F(cos 0, 0* )= 1, we have r(0 + )= (ro 3 + l) ,/: \ so the first thing to happen is that 
the fluid pops out about the origin to accommodate the injected fluid. Subsequently, because F(cos 
0,t)= 1 4- 0(t' 2 ) for small positive t, the initial displacements begin to change very slowly. 

We noted earlier that the velocity u is singular at x — but does not have a net in- or outflow 
there after t = 0. As a result, fluid is taken into the origin from some directions and given out in 
others. This is evidenced in the consequence of (4.5) that on the x$ axis all the points that satisfy 
r ^ 1 before the flow starts are drawn into the origin during the interval 7r/2 ^ t ^ it. Further 
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out along the * 3 axis displacements oscillate up and down without ever reaching the origin. It also 
follows from (4.5) that the displacements in the plane x 3 = grow without bound. 

Vertical displacements in a stratified fluid are of considerable interest, for they are often 
the easiest to measure. We note that they can be obtained equally well by interpreting (4.3). If we 
designate a particle in terms of its original height zo= r — cos 6 and #, the relation (4.3) also gives 
z= rcos 6 by 



2(20, 0; = 



zHr^cos 3 02F(cos<M) 
47rpo 



1/3 

(4.6) 



In a subsequent article [4] we shall give a more detailed analysis of the displacement and 
velocity fields of the fundamental and other singular solutions. 
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